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In the pursuit of speculative new particles, forces, and dimensions with vanishingly small influence 
on normal matter, understanding the ultimate physical limits of experimental sensitivity is essential. 

Here, I show that quantum decoherence offers a window into otherwise inaccessible realms. There is 
a standard quantum limit for diffusion that restricts some entanglement-generating phenomena, like 
soft collisions with new particle species, from having appreciable classical influence on normal matter. 

Such phenomena are elassieally undeteetahle but can be revealed by the anomalous decoherence they 
induce on non-classical superpositions with long-range coherence in phase space. This gives strong, 
novel motivation for the construction of matter interferometers and other experimental sources of 
large superpositions, which recently have seen rapid progress. Decoherence is always at least second 
order in the coupling strength, so such searches are best suited for soft interactions (e.g., small 
momentum transfers), but not weak ones (i.e., small coupling constants). 


Ultrasensitive detectors of feeble signals confront fun¬ 
damental quantum barriers that can only be beaten us¬ 
ing non-classical states of matter and radiation BEI- 
One way to capture the sense in which devices exploit¬ 
ing quantum superpositions can be strictly more sensitive 
than those using classical techniques is to derive a stan¬ 
dard quantum limit [I1[3H5] (SQL). In general, an SQL is 
a sensitivity limit arising from the assumption that the 
experimental probe used to investigate some phenomena 
can only be prepared and measured in certain ways — 
usually corresponding to the position basis — rather than 
the vast menu of quantum states and Hermitian observ¬ 
ables. For instance, the most common SQL restricts the 
detection of weak forces through sequential position mea¬ 
surements of a test mass, a limit that can be beaten by 
measuring a relative phase shift induced on a coherent 



FIG. 1. Decoherence detection with a Mach-Zehnder 
interferometer. A probe initially in a superposition jA'i) + 
|A^ 2 ) of two wavepackets is effectively measured by an interfer¬ 
ometer in the basis {|A^±) = lA/”!) ± |A^ 2 )}. If the probe scatters 
a hypothetical particle IF 0 ) into distinct conditional states 
\Ei) and IF 2 ) while being negligibly deflected, then the paths 
decohere and the particle will be found in state |A^-) with 
probability P- = (1 - Re {Ei\E2))/2. 


spatial superposition of the mass. This limit now chal¬ 
lenges state-of-the-art gravitational wave detectors like 
LIGO [2]. 

In this article, I derive an analogous SQL for detecting 
hypothetical weak sources of momentum diffusion, such 
as a soft flux of a new particle species ISHS], and show 
that this limit can be beaten using coherent superposi¬ 
tions that decohere rather than acquire a phase shift. I 
point out that such decoherence only enters at second or¬ 
der in the coupling between the source and the probe, in 
contrast to the measurement of unitary effects (like gravi¬ 
tational waves) which may be first order. Finally, I give a 
precise sense in which hypothetical phenomena, including 
both forces and diffusion, can be classically undetectable 
in an A ^ 0 limit while their influence on quantum ex¬ 
periments remains finite. 

Understanding the detective capabilities of experi¬ 
ments producing non-classical superpositions is a timely 
topic; bigger superpositions (as measured by mass, spa¬ 
tial separation, and lifetime) are more sensitive to de¬ 
coherence, and there has been spectacular recent suc¬ 
cesses in generating and verifying such states cni. Large 
molecules containing over ten thousand nucleons have 
been interfered in diffractive-slit experiments nn, with 
orders of magnitude of improvement predicted in the 
medium-term future cni E US]. Optically trapped 
nanoparticles promise to achieve wide spatial superpo¬ 
sitions of even greater amounts of matter [H [IS]— 
perhaps reaching an astounding 10^^ amu in spaceborne 
experiments mini —and quantum micromechanical res¬ 
onators will involve over 10^^ amu, albeit with femtome- 
ter superposition separations nans]- Such massive su¬ 
perpositions are predicated on the careful suppression of 
conventional sources of decoherence, making them un¬ 
usually sensitive to anomalous decoherence from novel 
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sources. General strategies exist to distinguish novel 
sources of decoherence from merely misunderstood con¬ 
ventional ones [9]. 

Framework. —Consider an interferometer (Fig. pro¬ 
ducing a coherent superposition |Ari) + |Ar 2 ) of some probe 
(e.g., a test mass), where |A^i) and lA/^ 2 ) are wavepackets 
separated by a large distance L. The wavepackets are 
recombined and, by way of appropriately aligned mirrors 
and detectors, the probe is effectively measured in the 
basis {\N^) = jA^i) ± |Ar 2 )}- Up to experimental error, the 
result is always |A^+). Now imagine that during its flight 
through the interferometer the probe interacts with some 
other hypothetical quantum system in an initial state 
\E 0 ) such that the hypothetical system is disturbed but 
there is negligible influence on either probe wavepacket 
individually. That is, assume that the evolution is well 
approximated by the form 

[ |iVi> + |iV2) ] lE^) |iVi) |i?i) + |iV2> \E2) (1) 

where \Ei) and \E 2 ) are arbitrary conditional states of 
the hypothetical system. 

If \Ei) and \E 2 ) differ only by a phase, {Ei\E 2 ) = 
then the probe remains in a pure state and, were it ob¬ 
served before the wavepackets are recombined at the end 
of the interferometer, would still be found with equal like¬ 
lihood in either arm. Nevertheless, the relative phase be¬ 
tween the two wavepackets can be inferred from the fact 
that the effective measurement at the end of the interfer¬ 
ometer will now result in outcome |A^_) with probability 
P_ = (1 - cos^)/2. This is in fact the basic mechanism 
by which atom interferometers are proposed to detect 
very weak forces like gravitational waves EDI, which are 
too small to displace atoms by a distance comparable to 
their wavepacket size yet still leave a measurable relative 
phase. 

But it’s also possible to consider scenarios for which 
(Pi 1 ^ 2 ) ^ 0, such as flux of soft particles scattering into 
distinguishable out states ISHS]. In this case, the probe 
is decohered and the probabilities for both outcomes of 
the measurement are equal: P+ = P- = 1/2. Just as for 
relative phases induced by weak forces, this decoherence 
can be detected even though the classical effects on the 
probe wavepackets (spatial displacements or momentum 
transfers) are completely negligible. 

Although searches for anomalous decoherence arising 
from speculative theories of quantum gravity m and 
modifications of quantum mechanics m have been con¬ 
sidered before, it is only recently that this technique has 
been proposed as a way to detect new particles within 
an otherwise conventional framework. In particular, cer¬ 
tain classes of low-mass (sub-MeV) dark matter — which 
would be invisible to traditional direct-detection experi¬ 
ments — could be identified by the characteristic deco¬ 
herence they induce in future matter interferometers [9]. 


Standard quantum limit .—Let us recall the popular 
heuristic argument for the SQL for the detection of 
weak forces [3]. A single quantum degree of freedom 
X, which we will imagine to be the position of a test 
mass m along one dimension, evolves in the presence 
of a small, spatially uniform force E for some time in¬ 
terval T. At t = 0, the mass is prepared by a posi¬ 
tion measurement of finite precision in a spatially lo¬ 
calized pure state with some finite, adjustable width 
^prep^ By the uncertainty principle, this prepared state 
necessarily has momentum width at least as large as 
= /i/(2crP^®P). During the time period T the mo¬ 
mentum width produces an additional position spread 
= f dTcrP^®P/m = hT/(2mcrP^®P) through wavefunc- 
tion dispersion. The total uncertainty for a final position 
measurement is then obtained by adding in quadrature: 
(^meas )2 _ j'^prep ^2 j'^disp^ 2 ^ ^ force translates the 

wavepacket uniformly [26] by Axp = FT^/(2m) relative 
to the null (F = 0) dynamics, and so the force will not 
be detectable if this distance is much smaller than 
One cannot simply make the initial position uncertainty 
arbitrarily small to increase sensitivity because of 
the resulting larger wavepacket dispersion ct^^^p. Mini¬ 
mizing with respect to the choice of initial width 

(jf^P yields o-”®"-" = A<tP’'®p = ^hT/m. We then obtain 
the requirement that a force must satisfy E > Fsql to 



FIG. 2. Standard quantum limit for forces and mo¬ 
mentum diffusion. A test mass is initially placed in a 
minimal uncertainty wavepacket with a Wigner distribution 
W{x,p) over phase space (top) that contains the bulk of its 
mass within a 2cr-contour of a Gaussian distribution (dashed 
black line). After a time T, the state will spread in posi¬ 
tion due to wavefunction dispersion even in the absence of a 
force or diffusion (solid black). The probability distribution 
of a subsequent position measurement (bottom) is obtained 
by integrating the Wigner distribution over momentum, (a) 
If a uniform force F displaces the wavepacket, it will not 
be reliably distinguishable from the zero-force scenario when 
F < Fsql (blue), in contrast to the case F » Fsql (red), (b) 
The standard quantum limit for momentum diffusion, Dsql, 
provides a similar threshold for the detection of an environ¬ 
ment that weakly smears the quantum state in phase space. 
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FIG. 3. Diffusion sensitivity from fine phase-space oscillations, (a) The Wigner distribution over phase space for a 
coherent superposition of two minimal uncertainty wavepackets for a particle free to move in one dimension. The wavepackets 
have position spread = yJhTJm^ average position xq = ±6crg^®^ (so L = 2xq = 12crJ’^®^), and average momentum po = 0. The 
grid units are and The fine (snh-h) phase-space oscillations near the origin are indicative of coherence 

m- (These oscillations are absent from an incoherent mixture of wavepackets, for which the Wigner distribution is given by 
the simple sum of the wavepackets.) (b) Under free evolution, the position and momentum of the particle become positively 
correlated, skewing the distribution. The additional effects of momentum diffusion over the time period can then be obtained by 
convolving this freely-evolved Wigner distribution with a Gaussian kernel ElES]; the 5cr-contour of this kernel for one choice 
of diffusion strength below the SQL is given by the black oval. This kernel is large compared to the oscillations, but small 
compared to the wavepacket widths, (c) The Wigner distribution after the initial state is evolved in the presence of momentum 
diffusion. The smoothing induced by the convolution only negligibly distorts the wavepackets, but the fine oscillations are 
destroyed, indicating that the state is almost fully decohered. Without the quantum coherence there would be no oscillations 
and hence very little sensitivity to the diffusion. The wavenumber of the oscillations is proportional to the separation L, so 
arbitrarily weak diffusion over a fixed time can be detected for sufficiently large separation. 


be detectable, where Fsql = 2yJhmlT^ is the SQL for 
weak-force detection with mass m acted upon for time 
interval T. See Fig. [^(a). 

Now we generalize to an open quantum system, allow¬ 
ing the mass to interact with an environment so that its 
reduced dynamics are generally not unitary. We seek to 
understand the physical limits for detecting the effects of 
an environment that weakly decoheres the test mass. For 
ideal quantum Brownian motion (QBM) in the friction¬ 
less limit, the density matrix of the test mass is governed 
by the master equation [iElEZlEB] 

dtPt = ^[x,[x,pt\\ (2) 

where H = j2m and where D is the momentum dif¬ 
fusion coefficient. One can get some intuition for the 
resulting decoherence by noting 0122 EH that the off- 
diagonal terms of the position-basis density matrix are 
exponentially suppressed for sufficiently large separation 
X - x'\ x') po{x, x') exp[-Dt(x - 

Equation is the simplest possible model of spatial 
decoherence [29], playing the pedagogical role of a har¬ 
monic oscillator for the study of open quantum systems. 
It is widely applicable due to the ubiquity of spatially lo¬ 
cal interactions, providing a good approximation to the 


effects of a thermal bath of linearly coupled oscillators 
[27] (although note limitations [30] |3T]). Importantly, 
(§ well describes the dynamics taken by a test mass 
subjected to coUisional decoherence 0 El ESI El from 
an environment of lighter particles, such as a gas 1331 , 
blackbody radiation [6] , or low-mass dark matter [9] . For 
a given environment, the approximation becomes more 
accurate for larger test masses. 

Intuitively, diffusion D that is small compared to other 
parameters will b e negligible, so like the force SQL 
(^SQL ~ yJhmjT^) we can guess the diffusion SQL up 
to a numerical factor by dimensional analysis: Dsql ~ 
This can be fleshed out with a heuristic ar¬ 
gument analogous to the one above [34]. Let the ini¬ 
tial state of the test mass be a wavepacket with spa¬ 
tial width = yJhT l(2m) minimizing (in the ab¬ 
sence of diffusion) The momentum diffusion in¬ 

duced by causes the test mass to follow a random 
walk in momentum space: - ‘^DT, yielding 

= f dT /m = VSDT^l{3m). The diffusion can 

only be reliably detected when There¬ 

fore, the diffusion SQL is that the diffusion strength obey 
D > Dsql for Dsql = 9hml{ST‘^). When D « Dsql, the 
diffusion cannot be detected. See Fig. [^(b), and the Ap- 
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pendix for a rigorous derivation. 

Like the force SQL, the diffusion SQL can be beaten by 
placing the test mass in a coherent superposition of two 
wavepackets widely separated in phase space (Fig.[^. To 
see this, consider an initial wavepacket |A^i) that min¬ 
imizes If D « Dsql then such a wavepacket 

is left essentially undisturbed by the diffusion over the 
timescale T. However, if a coherent superposition of 
|A^i) with |A^ 2 ) is prepared, where the latter is sim¬ 
ply the former translated by a sufficiently long distance 
L » hl\/DT^ then the superposition decoheres into an 
incoherent mixture of the two wavepackets within a de¬ 
coherence time td = h^l{DL‘^) « T. The wideness of 
the superposition, rather than its mere existence, is im¬ 
portant. The smallest diffusion coefficient that effec¬ 
tively decoheres a superposition of extent L is about 
Dmin = which reduces to the classical sensi¬ 

tivity Dsql when the superposed wavepackets are min¬ 
imally separated, i.e., when the separation L is of order 
the optimal wavepacket size = yJhT/m. 

No entanglement at first order .—It is worth noting 
that a probe Af cannot become entangled with another 
system S at first order in the Hamiltonian that cou¬ 
ples them, suggesting that anomalous decoherence is a 
method best suited for detecting interactions that are 
very soft (e.g., small momentum transfers) but not weak 
(i.e., small coupling constant). 

Let J\f and S initially be in an uncorrelated pure state, 
|'0o) = l^o) l-Eo), and suppose they evolve under Ut = 
ex-p(-iHt) where the Hamiltonian is of the form H = 
Hj\f + Hs + eHj. Here, ® Ig and Hg = <S) Hg 

are the self-Hamiltonians for J\f and S, respectively, and 
eHj is the interaction. Local unitaries cannot change 
the entanglement, so without loss of generality we can 
consider the modification U[ = ^+%Hst^-%Ht^ which 

peels off the unimportant local evolution of J\f and S. We 
will use the Zassenhaus formula [35] 

oo 

gA+S ^ gAgS J-j- gCnCA.S)^ 

n=2 

where Cn{A,B) is a homogeneous Lie polynomial of 
degree n with, for example, C 2 {A,B) = and 

C 3 {A,B) = l[B,[A,B]] + l[A,[A,B]]. Then, using 
A = -it{Hj\f + Hg)^ B = -iteHj, and the fact that Hj^ 
and Hg commute, we get 

U' = e-^^^‘*[l + 0{€^)] ( 4 ) 

for some new Hg that is independent of e, because each 
Cn{-it{Hj^ + Hg), -iteHj) is at least first order in e. So 
the probe’s density matrix is likewise pure to this accu¬ 
racy: 

p^rit) = Tvs [Ui\i>o){MUi^] = ly )(y I + O(e') (5) 

where |iVt) = {l - iet{Eo\Hi\Eo)) |iVo>. 


Decoherence of an initially pure state necessarily re¬ 
quires entanglement with an environment, so a signal 
of anomalous decoherence, which is linear in the off- 
diagonal elements of p^f{t), would be at least second 
order in the coupling. Therefore, interactions with 
very tiny coupling constants (e.g., the Fermi constant 
in searches for relic neutrinos) are poor candidates for 
detection through decoherence. Rather, an anomalous 
decoherence search is a method best suited for interac¬ 
tions with large coupling but with classical effects (e.g., 
momentum transfers or spatial displacements) that are 
negligible. Contrast this with traditional quantum en¬ 
hanced measurements of phase shifts, such as from weak 
gravitational waves, which are first order in the coupling 
strength (although they may still be small relative to de¬ 
coherence for other reasons ISI)- 

Classieal undeteetahility. —The exceptional sensitivity 
of superpositions to both displacement and decoherence 
can be attributed [23l|36] to fine structure in a test mass’s 
Wigner distribution at scales smaller than h. Indeed, 
the decoherence of a superposition of widely separated 
wavepackets is the “complexification” of the phase shift 
induced by a force. Let us consider the class of completely 
positive (CP) reduced dynamics for the probe’s den¬ 
sity matrix that leave two wavepackets undisturbed, as 
in Fig. for j = 1,2. One can 

check that any CP map satisfying this requirement must 
operate on the two-dimensional subspace spanned by the 
(essentially orthogonal) wavepacket states as 


4>[p] = 4> 


/ Pii 
\ P21 


Pl2 \ ^ / Pii 
P22 ) \ 7*P21 


1P12 \ 
P22 /’ 


(6) 


where j = {Ei\E 2 ) = e is the deeoherenee faetor, with 
0 real and s > 0. 

The classical force obeys 0 = ELTjh and (since it 
is unitary) s = 0. Decoherence from momentum diffu¬ 
sion, (§. obeys s = DLfTjhf and (since it is isotropic) 
^ = 0. These can be obtained while still satisfying 
Ri |A(^ )(A(^ | to arbitrary accuracy by holding 
L, T, Ejh, and Djhf fixed and taking ^ 0, so that 
E,D ^ 0. In other words, in the classical limit where 
wavepackets become points in phase space |3Zl|38] (i.e., 
states of arbitrarily precise position and momentum), the 
phase and decoherence induced by $ in an interferometer 
can remain fixed — and hence detectable — even as the 
classical strength of the force E and the diffusion D go to 
zero. Similar limiting behavior exists if the force or dif¬ 
fusion have more complicated dependence on position, or 
if both effects occur simultaneously (e.g., the anisotropic 
dark matter “wind” [9]). It is in this sense that coherent 
superpositions enable the detection of phenomena that 
are elassieally undeteetahle. 
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Ideal QBM and contractive states 

It is well known that the popular heuristic argument 
for the standard quantum limit (SQL) for weak forces 
[n [sHS] contains a loophole [39lfi3] . and this issue ap¬ 
plies likewise to the argument for the diffusion SQL. 
Even when restricting to the simple case of a uniform 
force acting on an initial Gaussian state of the probe 
that is sensibly described by its first and second mo¬ 
ments, a counter-example is obtained by preparing the 
probe in a “contractive” (but still Gaussian) state for 
which position and momentum are negatively correlated 
[39l [43]. (Contractive states can equivalently be defined 
as squeezed coherent states such that the line in phase 
space along which the state is positively squeezed has 
negative slope.) Like coherent superpositions, contrac¬ 
tive states are non-classical and very difficult to prepare 
experimentally. A rigorous SQL can be obtained in the 
general case of ideal quantum Brownian motion (QBM), 
which encompasses both weak forces and diffusion, when 
non-contracting Gaussian initial states are assumed. 

First let us place our model of momentum diffusion and 
spatial decoherence in a general framework of ideal QBM. 
The simplest type of open quantum system dynamics are 
Markovian and time-homogeneous, forming a quantum 
dynamical semigroup for which the purity of the state 
is nonincreasing with time mms]. Setting h = such 
evolution is described in generality by a Lindblad master 
equation 

dtpt = -i[H,pt\ +y . 

(7) 

with the unitary part generated by the Hamiltonian H 
and the non-unitary part by the Lindblad operators 
The smallest class of Hamiltonians containing the har¬ 
monic oscillator and invariant under rotations and trans¬ 
lations in phase space are those that are polynomials no 
more than quadratic in the phase space variables x and p. 
Gonstant terms in Lindblad operators can be absorbed 
into the Hamiltonian so, to generalize to open quantum 
systems, we consider the simplest nontrivial Lindblad op¬ 
erators: those linear in x and p. Linear terms in the 
Hamiltonian can then be eliminated by translating coor¬ 
dinates X ^ X + xq, p ^ p + po, putting the origin at the 
center of the ellipse or hyperbola defined by H. (This 
can’t always be done in special cases when quadratic 
parts of H are zero, such as for a uniform force, which 
must be handled separately EHO Using Einstein summa¬ 
tion (a, 6 ,... = x,p) over the two phase-space coordinates 


and letting a = (x,p), we parameterize the dynamical op¬ 
erators as H = 1^^ and where Hat is 

a real symmetric matrix and the are complex vectors. 
We can then change variables to Dab = Re )* 

A = = ImE,(LP)-^LP, and Kab = 

Hab + where is the antisymmetric Levi-Givita 
symbol with = +1 = 621 - This specializes 0 to QBM: 

dipt = [a“, {a^ pt}] - Ua 6 [«“, [a^ Pt]] ■ ( 8 ) 


We move to the Wigner representation, where the 
quantum state pt is exchanged for a quasiprobability dis¬ 
tribution Wt{x,p) over phase space: 


Wi 


(x,p) = j 


dAx 


-ipAx 


, Ax Ax^ 
Pt{x+—,X-—) (9) 


such that fdxWt(x,p) = ptip,p) and fdpWt{x,p) = 
Pt{x,x). Under 0. the Wigner distribution correspond¬ 
ing to the state pt obeys a Klein-Kramers equation 

[MEH] 


dwt{a) = [-K\daa'^ + ]^D‘''’dadb]Wt{a) (10) 

where a = {x,p), da = djda^^ and = 

The only terms that generate non-Hamiltonian evo¬ 
lution are the real parameter A = and the real 

symmetric matrix . The matrix must [28] have 
non-negative eigenvalues and satisfy \D^^\ > A^, where 
= DabD^^l2 is the matrix determinant. One can 
always diagonalize through a coordinate rotation in 
phase space, and the resulting diagonal elements quan¬ 
tify the diffusion in the new x and p directions; by the 


symplectic symmetry of ( 10 ), they have mathematically 
similar behavior. 

Although it is possible to consider general SQLs for 
measuring all these parameters, they are closely related 
by symmetries and our primary concern is with the spe¬ 
cial case of pure momentum diffusion discussed in the 
main text (|^. This is when = 1/m, = 2D, and 

all other coefficients of and are zero. In the 


Wigner representation, (10) reduces to 




m dx dp^ 


Wt(x,p), 


( 11 ) 


which is identical in form to the Klein-Kramers equa¬ 
tion for a classical phase-space probability distribution 
exhibiting diffusion in momentum space, i.e., ideal clas¬ 
sical frictionless Brownian motion [28]. (This is very dif¬ 
ferent from the large-friction limit studied by Einstein 
and Smoluchowski.) 

The exact solution to ( pT) ) for any initial Wigner dis¬ 
tribution Wo{a) is [Ml[25] 


Wt{a)=gt{a)^Wo{R-t{a)) 


( 12 ) 
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where R-t{a) = R-t{x,p) = (x - ptlm^p) is the reversed 
Hamiltonian evolution, ★ denotes the convolution opera¬ 
tor in phase space, and gt{<^) ^ exp(-(C'^“^)“^aaa 5 / 2 ) is 
the normalized Gaussian smoothing kernel with a time- 
dependent covariance matrix given by 


i/(2m) 1 ) 


(13) 


The distribution of outcomes for the final position mea¬ 
surement is obtained by integrating the Wigner distribu¬ 
tion over momentum: P(x) = f dpWrix^p). 

If contractive states are excluded, one can check that 
the Gaussian initial states producing final measurement 
distributions depending most sensitively on D (in the 
sense that the Ghernoff bound exponent is maximized 
[46]) are always the ones for which position and momen¬ 


tum are initially uncorrelated; in this case, the Heisen¬ 
berg uncertainty inequality is saturated and the uncer¬ 
tainties (cr’s) from the heuristic arguments in the main 
text coincide with the standard deviations of the respec¬ 
tive Gaussian probability distributions up to factors of 
order unity. A similar SQL for weak forces m can be 
proved with the same techniques using H = p^ j2m - Fx. 

Note again that K\ and were assumed to be time- 
independent, consistent with time-homogeneity and the 
Markov property, and well motivated by the observation 
that the time-dependent components of these parame¬ 
ters in real-world cases of quantum Brownian motion 
are often transient when sourced by thermal baths m- 
Gare must be taken with this assumption [30l [31], but 
more complicated SQLs can be derived even when it is 
relaxed. Likewise, one can straightforwardly derive tra¬ 
ditional SQLs for time-varying forces or effective masses. 


